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ABSTRACT 

An analytic solution to the problem of radiation from an open-ended 
parallel-plate waveguide into a dielectric or plasma sla£ is obtained. 

The electromagnetic fields are formulated in terms of a discrete mode 
spectrum inside the waveguide and continuous mode spectra in the open 
regions. Mode matching at the waveguide aperture plane leads to a sin- 
gular integral equation. A solution of this equation allows the deter- 
mination of the fields throughout the whole of space. However, this 
equation is not solved directly, but : instead, an auxiliary integral equa- 
tion is formed from the original equation by employing the known exact 
solution of the canonical problem, viz., radiation from an open-ended 
waveguide into free-space. The auxiliary equation has the advantage 
that it may be solved exactly in an asymptotic sense. This asymptotic 
solution demonstrates that the correct edge condition is satisfied. 
Standard numerical techniques are employed, and results giving the fields 
and distribution of scattered power for a TEM incident mode are presented. 

Solutions for other problems are also formulated, viz., the flanged 
waveguide and a guide radiating into a reflecting screen. These prob- 
lems may also be solved when the geometry ig circular. Furthermore, 
the edge condition is also shown to be satisfied in circular geometry. 
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1. INTRODUCTION 


Exact solutions can be found for certain open-region boundary value 
problems in which the geometry can be separated into several distinct 
semi-infinite regions. For example, an open-ended parallel-plate waveguide 
radiating into a homogeneous space has been solved by Noble (1958) using the 
Weiner-Hopf technique and by Mittra and Bates (1965) employing an extension 
of the function- theoretic technique to open-regions. Recently, problems 
in which one of the regions is modified have also been found solvable. 

Bates and Mittra (1967), using the Weiner-Hopf technique, have found an 
analytic solution for the problem where a dielectric or plasma slab is 
excited by a semi-infinite waveguide (cf» Figure 1). An important 
characteristic of this problem is that the medium filling the open-region 
is non-uniform in the transverse direction, while the complete geometry is 
longitudinally uniform. When the geometry is non-uniform in the longitudinal 
direction, a class of boundary value problems results for which the Weiner- 
Hopf technique is no longer directly applicable. Figure 2 shows a 
longitudinally non-uniform problem, i.e., a semi- infinite parallel-plate 
waveguide radiating through a dielectric or plasma slab. The purpose of 
this work is to present a method based on the function-theoretic technique 
which is useful in attacking problems which have longitudinal non- 
uniformities. 

Formulation of a problem by the function-theoretic technique is 
typically accomplished by employing the characteristic or normal modes 
in each region. Such a modal representation can be used to advantage 
when there is a longitudinal inhomogeneity within an individual region. . 
This applies to stratified media as well. Such regional non-uniformities 
can be taken into account by the application of conventional mode-matching 
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procedures at the interfaces between the various layers or subregions. 

An application of the continuity conditions at the aperture-plane of the 
waveguide, leads to an integral equation, which is an extended form of the 
corresponding equation for the unmodified geometry. An exact solution of 
this equation for the canonical problem, i.e., for the unmodified geometry, 
is available (Mittra and Bates, 1965) and is used in the construction 
of the solution of the problem shown in Figure 2. 

A function F(co) of a complex variable u> is constructed having a 
certain pole-zero structure and specified branch singularities. Certain 
manipulations with this function in the complex plane yield an integral 
equation which is then compared to the one obtained by mode-matching. By 
choosing a particular form for F(w), together with the above comparison, 
an auxiliary integral equation is formed. This auxiliary equation may be 
solved instead of the original integral equation. 

There are two definite advantages gained by employing the auxiliary 
integral equation . These advantages are: (1) The auxiliary equation is 

in a more efficient form for numerical solution than is the original 
integral equation. (2) An exact solution in the asymptotic sense of the 
auxiliary equation is possible. It is this asymptotic behavior which 
demonstrates that the correct edge condition is indeed satisfied 
(Meixner, 1954). It is often difficult if not impossible to prove the 
satisfaction of this condition when alternative methods are employed in 


constructing a solution 
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2. RADIATION FROM AN OPEN-ENDED PARALLEL-PLATE WAVEGUIDE INTO AN 

INHOMOGENEOUSLY FILLED SPACE 

2„1 Formulation of the Problem 

Consider an open-ended parallel-plate waveguide radiating into an 
open-region which is partitioned by a dielectric slab. The geometry of 
this structure is shown in Figure 2. The waveguide is excited from the 
left by a single TEM or TM mode (p = 1,2,....) of unit amplitude in 
region B. The space surrounding the guide is divided into four regions: 

A, C, D, and E. A right-hand Cartesian coordinate system is placed at 
the guide aperture with the pertinent dimensions given by b, the guide 
half-width; the guide-aperture separation; and t, the slab half-width. 
The relative permittivity of region A is given by k while that of the 
remaining space is taken as unity. 

A solution of the source-free Maxwell’s equations for this geometry 
leads to a discrete eigenmode expansion of the electromagnetic fields in 
region B and transform representations in terms of the continuous mode 
spectra in the remaining open-regions. The transverse fields are matched 
at each of the regional interfaces, and the resulting equations are 
manipulated in order to yield expressions for the mode and spectral 
weight coefficients. These coefficients are then explicitly solved for 
by an appropriate modification of the function-theoretic technique. This 
technique, in its original form, was first applied to closed waveguide 
problems by Hurd and Gruenberg (1954) and to open periodic structures 
by Berz (1951) and Whitehead (1951). These early works yielded infinite 
matrix equations, the elements of which were functions of the longitudinal 
wave numbers and mode coefficients for the various regions. Mittra and 
Bates (1965) extended the function-theoretic technique to an open-region 
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problem by employing a limit process to the associated closed-region 
problem. Due to the continuous nature of the eigenspectrum in the open- 
regions, an integral equation results corresponding to the infinite matrix 
equation for the associated closed-region problem. For the modified 
problem of Figure 2 s an extended version of this integral equation is 
obtained, which can be transformed into a matrix equation by the usual 
numerical techniques. 

The geometry of Figure 2, being y-independent and x-symmetric, 
scatters only TEM and TM type waves with respect to the x-£oordinate . 

The total magnetic intensity $ = in each region is represented as 
follows : 

$ = — 

A 2 

a 

l < z < £ + 2t, (2.1) 


A^(a) cosh[x (z-£-t) ] + A£(a) sinh[r (z-£-t) ] 


exp (-j ax) da s 


CO 

$ B = cosp”j exp (-8 z) + Z [B g cos|-™-j exp(3 g z)], 

“ - ” s=o ' 

fxj <_ b, z < 0, (2.2) 


$ 


C 



C(y) exp(nz) exp [ - j y (x-b ) ] dy, 


x > b, z <_ 0, 


(2.3a) 


*C 



C(y) exp(nz) exp[-jy (x+b) ] dy. 


x <-b, z 0, 


(2.3b) 
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% = J [jD + (a) exp(-£z) + D (a) exp(£z)J exp(-jax) da, 
a 

0 < z < £, (2.4) 

and 

$ E = -| | E(a) exp [-g (z-£-2t)] exp (-j ax) da, 
a 

Jt + 2t £ z. (2.5) 

The mode coefficients are given by B , (s = 0,1,2...), in region B, and 

s 

the spectral weight functions by A^(a), A2(a), C(y), D + (a), D (a) and 
E(a) in their respective open-regions. The fields have a traveling wave 
representation in each region except A, where the representation is 
in terms of even and odd modes with respect to the z-coordinate . This 
is a convenient representation for region A when there are surface modes 
present within the slab. The variables of integration a and y lie on 
the path a shown in Figure 3. Use of the exp(jcot) time convention 
requires that the following function branches be employed in order to 
yield outgoing waves: 


a/ 2 2 

j/k -a 
J o 




(2.6b) 
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The free-space propagation constant is given by k Q = where u 

is the angular frequency and and e q are respectively the permeability 
and permittivity of free-space. The branch cuts for (2.5) through (2.8) 
are shown in Figure 3 for K-real and greater than unity. The remaining 
components of the electric field are obtained from the following 
relationships : 


_ v-1 3* 

E x - 0** o 'r> al 


( 2 . 9 ) 


and 


E 

z 


Owe e ) 
J o r 


3 $ 

9x 


( 2 . 10 ) 


where e^ = k in region A and unity throughout the remaining space. 

2.2 Representation in the Spectral Domain 

The first step toward obtaining the solution of (2.1) through (2.5) 
is to express the unknown mode and weight coefficients in terms of 
E(ot), the weight coefficient for the forward open-region. The procedure 
is to enforce the continuity Conditions of the transverse field components 
at the three regional interfaces. Continuity at z=0 gives 


$ B (x,o) = $ D (x,o) 


( 2 . 11 ) 


H x l 1 b 


9$ b (x,o) 9$ d (x,o) 


dZ 


9z 


( 2 . 12 ) 



10 


and 


$ G (x,o) = $ D (x,o) 




> M 1 b 


34> c (x,o) 9 $ d (x,o) 

3z 3z 


Similarly, at z = £ we obtain 


$ d (x,£) = $ a (x,£) 

3$ d (x,£) 3$ a (x,£) 

3z K 3z 

while at z = £ + 2t we have 


$ A (x,£+2t) = $ E (x,£+2t) 

and 

_ 1 3$ A (x,£+2t) 3$ E (x,£+2t) 

K 3z ~ = 3z ° 

Combining (2.17) and (2.18) yields the expressions for the even 
odd spectral weight functions in region A 


A^(a) = E(a) [cosh(tt) + ” sinh(tt)] 

and 

A£ (a) = E(a) [sinh(-rt) + cosh(xt)]. 


{2. 13) 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 

and 


(2.19) 


( 2 . 20 ) 
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Similarly, a combination of (2.15), (2.16), (2.19), and (2.20) gives 
the forward and backward traveling wave spectral weight functions in 


region D 




D + (a) = E(ct) R* (a) 

(2.21) 

and 




D~(a) = E(a) Q’ (a) , 

(2.22) 


where R' (a) = R(a) exp(5£) and Q’(a) = Q(a) exp(-££), with 


R(a) = cosh(2xt) + y 


El + _L 

T KE, 


sinh(2xt) 


(2.23) 


and 


Q(a) 


1 

2 


El 

X 


X 

<5 


] 


sinh(2xt) . 


(2.24) 


The above procedure, used to obtain expressions for the spectral 
weight functions in the regions lying between the guide aperture and the 
forward open-region, is not limited only to problems having two 
intermediate regions. One can express the weight functions in each 
region of an n-layered medium in terms of the corresponding function in 
the far-forward region. This procedure is similar to that of determining 
the voltage along each segment of a cascade of transmission lines, each 
portion having a different characteristic impedance, in terms of the 
voltage at the line’s termination. 

Fourier analyzing (2.11) and (2.12), i.e., multiplying each side 
of these relationships by cos(irnrx/b), (m=0,l, . . . . ) , and integrating with 
respect to x from -b to b, together with the relationships (2.21) and 
(2.22), yields 
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(-1) P b (1+6°) 6 P + (-l) m b B (1+6°) = 
pm mm 


E f (a) 

1 

3 

S— ' 

O' 

+ 

3 

P4 
i 

,2 .2 

a 

l_ m — > 


(2.25) 


and 


(-1) P b 8 (1+6°) 6 P - (-lfb 8 B (1+6°) = 

p v p m m m v m 


C E' (a) 


R* (a) - Q* (a) 


2 „2 

r-r 

m 


da. 


(2.26) 


where E* (a) = a sin(ab) E(a) and 6^ is the Kronecker delta defined by 




°> if i I s j 
1, if i = j 


(2.27) 


The integrands of (2.25) and (2.26) have poles at £ = + B ffi or 
equivalently when a = _+ imr/b. These poles are cancelled by the 
corresponding zeros of sin(ab) , which is a factor of E' (a) . Due to this 
pole-zero cancelation, the integrals of (2.25) and (2.26) need not be 
considered in the Cauchy principal value sense. 

Multiplying (2.25) by 8 m and first adding and then subtracting it 
from (2.26) yields 


(-l) p b 8(1+6°) 6j? 

p pm 


E* (a) 


R’(cO _ Q»(a) 


C-B. 


m 


S+B 


da 


m 


(2.28) 


and 


- ("D m b 8 B (1+6°) 
mm m 


E'(a) 


R* (a) Q* (a) 


?+B. 


m 


5-B, 


da 


m 


(2.29) 
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Note that the path a employed in all integrations previous to (2,28) 
has been reduced to o' in (2.28) and (2.29). This new path is identical 
to that portion of a which extends from zero to infinity and lies entirely 
in the upper-half a-plane (cf. Figure 3). Such a reduction is carried 
out by requiring E(a) to be a symmetric function of a, the remaining 
part of the integrand being obviously symmetric. 

Multiplying (2.13) by exp[jy(x-b)] and integrating with respect to 
x from b to «■, and repeating the same process with a multiplier 
exp[ jy (x+b) ] and integrating from -«? to-b yields 


it C (y) = it E(y) [R* (y) + Q'(y)] cos(yb) 


and 


P.V. 


E'(a) 


R* (a) + Q' (a) 

2 2 
? -q 


da 


IT n C (y) = - If n E(y) [R'(y) - Q'(y)] cos(yb) 


(2.30) 


+ P.V. 


E'(a)£ 


R'(a) ~ Q’(a) 

_2 2 
S -q 


da 


(2.31) 


Due to the continuous nature of q, a cancelling of the poles at 
E, = + q is not found except when q =,+ 8 . Therefore, the integrals 
in (2,30) and (2.31) are written as the Cauchy principal value type. 

Combining (2.30) and (2,31), as was done in obtaining (2.28) and 
(2.29), gives 


it q E(y)R ? (y) cos(yb) = P.V. 

» 

E'(a) 

R* (a) O’ (a) 

5-q ~ C+q 

a' 

. mm 


it q E(y)R ? (y) cos(yb) = P.V 


da (2.32) 
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and 


ir n C(y) tt n E(y) Q'(y) cos(yb) = 


P.V. 


E’(a) 


R* (a) _ O' (a) 

5+n C-n 


da . 


(2.33) 


Equations (2.29) and (2.33) express the two remaining unknowns, viz., 

B and C(y), in terms of E(a). Comparing (2.28) with (2.32), i.e., 
s 

when n = 3 , we obtain the values of E(a) at a countable-infinite 
m 

number of points 


E(sTr/b) = \ 


b (1+5°) JjitR 1 (pir/b)j ; 




(2.34) 


; s ^ p 


Equation (2.32) together with (2.34) represents a homogeneous integral 
equation for E(a) with prescribed function values at particular values 
of a. In the following section a method is presented to solve this 
integral equation based on an extension of the function-theoretic technique 
for open-regions. 


2.3 Modification of the Function-Theoretic Technique 

As stated in the introduction, a closed-form solution to the 
problem associated with the geometry of Figure 2, but with the dielectric 
or plasma slab removed (hereafter called the canonical problem) , has been 
obtained. When the slab is inserted, one finds that by employing a 
suitably modified version of the function-theoretic technique, an analytic 
solution is still possible. However, the solution is no longer in a 
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closed-form, but is most easily obtained from the solution of an auxiliary 
integral equation. 

2.3.1 Solution in the Spectral Domain 

Consider a function F(w) of a complex variable a), which is 
factored as follows : 

F(ui) = F 1 (oj) T(w) . (2.35) 

The factor F^(w) is that function which one employs in the solution of 
the canonical problem (Mittra and Bates, 1965) by the function- theoretic 
technique, and is given by 


F 1 (m) = K n(u,6) N(u>) (u-3 ) 


(2.36) 


The constant K is a normalization factor which depends upon the 
amplitude of the incident waveguide mode. 


n(o>,3) = n 

s=l 


; 

'i-f-l 


3 


si 


exp(a)/ 0 ) 


(2.37) 


is an absolutely convergent infinite product expansion, and 


N(w) = exp 


— [1-C + ln| 

tt e 1 


2tt 

k b 
o , 


i • bio 
] " 3 T + 


- /) 2 +k 2 In 

IT O 


1 ,/2]T 2> 

o>— /to +k 
c 

-j k 


(2.38) 


is obtained from the limit function of the ratio of two infinite product 
expansions associated with regions C and D in the canonical closed-region 
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problem as the outer guide is made infinitely wide (cf. Section 2.3.2 ). 

The constant C g in (2.38) is equal to Euler's constant. 

For large values of w, the, canonical function F^(ui) decays 

exponentially in the upper half w-plane and decreases algebraically as 

“1 / 2 

| to | in the , lower half-plane. Furthermore, F^(to) has a branch point 
at m = g Q and a branch singularity, shown in Figure 4, due to N(u>). 

As will become apparent in the sequel, the composite function F(ijo) must 
have a simple pole at to = g^, the propagation constant for the incident 
guide mode, and a branch point at m = - g . Therefore, we require 
T (to) to possess these added singularities with a branch cut along the 
reflection of the cut for F^(u>) (cf. Figure 4. ). The edge condition 
(Meixner, 1954) requires that F(cd) -+O0* 5 V ) us | oj | -*■ °°, arg(u)) 

= - it, and therefore, it is stipulated that T(w) V )> where 



The branch of arcsin is chosen as follows: 


(2.39) 


— — arcsin(0) <_ -jp j -1 <_ 0 <_1 . (2.40) 


The integral equation (2.32) for the weight function E(a) can now 
be, solved. Integration of the function F(co)/(w-n) along the path E 
(cf. Figure 4), with n placed directly on the branch singularity, yields 
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The integrals along the contours a and are zero in (2.41) due to 
the exponential decay and algebraic decrease respectively of the integrands 
along these paths as |co| °°. The integrals along the paths and a ^ 
are interpreted as being of the Cauchy principal value type. Employing 
the relationship 


N (to | ) = N (oj I ) exp (- 2jb t4 2 +k^) , 


(2.42) 


which expresses the discontinuity in N(to) across the branch cut, and 
by transforming the integral along a ^ and a ^ into a single integral 
along via the transformation of variable <a -*• -to, (2.41) is reduced to 
the following form: 


F(n 



exp(-jyb) 


cos(yb) = P.V. 



F (-to 


) - F(-o> l ) 
°3 ' 4 


2trj (w+n) 


dto . 


(2.42) 


The integral in (2.42) is next transformed into an integral along 
the path o', which was employed in (2.32), via the transformation of 
variable 01 = C and dw = a da/S* The result is given by 


F(n I ) exp(-jyb) cos(yb) = P.V. 
a l 


'a F(?| ) sin(ab) exp(-jab) 

g l 

5 ir (S-n) 


+ 


a [F(-?| >-F(-?| )]“ 


5 2-rrj a+n) 


(2.43) 
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A direct cpmparison of the integrands of (2.43) arid (2.32) yields 

F(c| a ) exp (-jab) 

E(a) ■ ■ — R -- fer» - ; — «•“> 


and 


E(a) 


[F(-?| ) - F(-£| o )] 

_3 a 4 

Q'(a) 2 tf j sin(ab) g 


A comparison of the non- integrated terms yields 


(2.45) 


F(n| ) cos(yb) exp(-jby) = it n E(y) R* (y) cos(yb) . (2.46) 

1 

Equations (2.44) and (2.46) are identical when n = g and y = a. This 
identity shows the self-consistency of the solution for E(a) as 
determined by (2.44) or (2.46). Furthermore, E(a), as obtained from 
(2.44), is seen to satisfy . the conditions at a = (sn/b), given by (2.34). 
It is from the value of E(a) at a *= (pir/b) that the normalization 
constant K employed in (2.36) is determined. With the aid of (2.34) 

and (2.38), K is then obtained from the relationship 

\ 

F(el )-b«„ <-l) p (1+6°). (2.47) 

pl o 1 p p 

In order that (2.44) and (2.45) be consistent expressions for E(a) 
they are first equated to one another, then the factored form of F(o>) 
is substituted, and finally the results are rearranged in order to 
obtain 

T(-gi ) “ T (“S L ) - 27rj X(0 T(C), 

°3 °4 


(2.48) 
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where A(£) is given by 


Q'(a) F-,(5| rt ) sin(ab) exp (-jab). 



(2.49) 


Equation (2.48) represents an expression of the discontinuity in T(-£) 
across the branch cut, as shown in Figure 4, in terms of the function 
value T(5) at the reflection point £. In Section 2.4, (2.48) is 
transformed into the integral equation which was previously referred to 
as the auxiliary integral equation. 

At this point, the problem is solved contingent upon the construction 
of T(co). However, simpler expressions exist for the mode coefficients 
B m and weight function C(y), than the integral representations given by 
(2.29) and (2.33). Integrating along the contour I with -8 m replacing 
n in (2.41), together with the subsequent transformation of variable 
and a direct comparison of the result with (2.29), yields 


B = 
m 


<-ir fhj . ) 

m o_ 


b (1+6°) 
m m 


(2.50) 


The same process with n replaced by -q yields 



Before Continuing on to the solution of T(u>), we digress so that 
we may compare the solution just derived with that for the corresponding 
closed-region problem. This is done to gain some insight into the form 
of the auxiliary integral equation. 
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2.3.2 Comparison with the Closed-Region Problem 

The closed-region geometry associated with the problem just solved, 
is shown in Figure 5. The outer guide now has a finite half-width given 
by a, In all of the regions, the fields are now represented by eigenmode 
expansions due to the discrete nature of the eigenvalue spectrum for 
closed-regions. The formulation of the problem follows closely that of 
Section 2.1, with the exception that all of the integrals are replaced 
by infinite summations, and the continuous propagation constants take 
on discrete values 


and 


4t) 2 - k o 

(2.52) 

4f) 2 - « < 

(2.53) 

At) 2 - < * 

(2.54) 


where c = a - b. 

Mode matching in the plane of the small guide's aperture and 
subsequent manipulations of the resulting equations, yields a.doubly- 
infinite set of matrix equations for the mode coefficients E g . This set 
of equations is given by 


(-1) P b g (1+6°) 6 P = 2b B E R* <5° + E E ? 

p p m ooom .,s 

r s=l 



s m 


% 


s +& 

s m 


(2.55a) 


and 


0 = 2c n E R 6 
ooom 


s=l 


E' 

s 


R* 

s 


5 - ir 

s m 


5 +n 

s m 


(2.55b) 
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Figure 5 


The closed region problem. 
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where 


R' 

s 




% 



and E* 
s 



E . 
s 


In the infinite summation of (2.55b), it is assumed that £ ^ n for all 

values of s and m. This is equivalent to requiring that the ratio of 
the outer and inner guide half-widths, viz,,a/b , is not a rational 
number. However, when this ratio is a rational number, we exclude from 
the summation in (2.55b) the q-th term when = n , and add to the term 
outside the summation the limit of the s = q-th term as s -> q» The 
appropriate limit is given by 


lim 
s -> q 


E' R' 
s s 




E R' n b 7 cos (y b) , 
q q q a q 


(2.56) 


where. we have written y = (qir/c). When the outer guide becomes 
infinitely large, i.e., c/a -»■ 1, a - c = b, (2.56) becomes the term 
outside the integral in (2.32) for the corresponding open region-problem, 
viz., - , rrTiE(Y)Rly) cos(yb). Due to the continuous nature of the eigen- 
spectrum in the open region case, this non- integrated term is always 
present, since we are forced to let £ = n and a = y. 

A solution of the doubly-infinite set of matrix equations is carried 
out by a modified version of the function- theoretic technique. The 
unmodified technique was employed by Hurd and Gruenberg (1954) for the , 
solution of matrix equations associated with the bifurcated waveguide 
problem. A form for modifying the function-theoretic technique so that 
it can be employed to solve the bifurcated waveguide problem, with one 
of the interior regions modified was first presented by Mittra, Lee, 
and Vanblaricum (1968). In this form, the infinite product expansion 
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associated with the modified region was written as a new infinite product 
expansion containing a countable-infinite number of shifts in the zeros 
of the original expansion, i.e., shifts from the zeros associated with the 
unmodified problem,, This zero-shift method is useful for closed region 
problems, but becomes awkward in the open-region case. The difficulty 
with the zero-shift formulation is that one cannot easily determine the 
limit function generated by the shifted infinite product expansions as the 
outer guide width becomes infinitely large. An alternative form for 
modifying the infinite product expansion is to employ an equivalent 
partial fraction expansion 



1 

/ 



°° 


1 “■ 

1 O 


= n(ui,n) 

1 1 , , y ® 


n +A 

exp n +A 


-L + 00 L 

, o-n 

_ 

s s 

Is sj 



S=1 s 


(2.57) 


Equation (2.57) shows the expansion containing the unknown shifts A on 

s 

the left, while an equivalent form is shown on the right. This equivalent 
form is written in terms of the unmodified expansion n(u>,n) and the 
partial fraction expansion coefficients g g . When the limit process is 
considered, i.e., (a,c) °=> and a - c = b, it is found that the partial 

fraction expansion has an integral form for a limit function 


1 + 0 


Z 

s=l 


o>— n 





1 + 0 


. 


a 


g(z) dz 
co-n(z) 


(2.58) 


The ratios of the infinite product expansion associated with regions 
C and D of. Figure 5 are found to have the following limit function: 


lim 

a,c -* 00 
a-c=b 


n(u),n) 

n(u>,5) 


-*■ N(co) . 


(2.59) 
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The function N(oj) is closely related to that given by (2.38). This 

relationship will be shown in Section 4. This function N(m) was 

originally obtained by the limiting process by Bates and Mittra (1965) . 

The solution of (2.55a) and (2„55b) for E by the function-theoretic 

s 

technique, is accomplished by the construction of a meromorphic function 
H(uj). This function in a factored form is written 


H(oj) = H^(w) 



Z 

s=l 



(2.60) 


where H^(oj) is that function employed in the solution of the canonical 
problem (Hurd & Gruenberg, 1954). Manipulations in the complex plane 
similar to those for the open-region problem (there is no branch 
singularity to circumvent in the contour integration) yield solutions 
for the various mode coefficients. 


B 

o 



2b e 0 


(2.61) 


(-l) s H(-B ) 

B g - - g-g — , (s = 1,2,...), (2.62) 


C 

o 


b 

a 



+ 


2c n o 


(2.63) 


C 

s 


W-n s ) 


, (s 1,2, 0 




(2.64) 


and 


E 

o 


b 1_ 

a R* 
o 


6 ° 

P 


(2.65) 
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residue H(co = £ ) 



(s — l,2,o..). 


( 2 . 66 ) 


E 

s 


residue H(u=-£ ) - 

s 1 



(s = 1,2, . . .) . 


(2.67) 


Elimination of E from (2.66) and (2„67) yields an infinite auxiliary 
s 

matrix equation for the expansion coefficients 


g = X 
6 s s 


1 + Z -!s_ 

=1 c +c 


Z -6 
s p 


m=l s m 


where 


(s - l,2,o..), (2.68) 


X 

s 


% 

R' 

s 


residue H, (w=£ ) 
1 s 

v-v 


(2.69) 


The normalization condition, which depends on the amplitude of the 
incident guide mode, is given by 


f 


H «p> ’ 


— b fl 
a o 




p=o 


y 


(-l) p b 6 , p^o 


V 


(2.70) 


As the outer guide becomes infinitely wide, we expect (2.68) to become 
the auxiliary integral equation associated with the open-region problem. 
This is indeed the case, and the discussion in the following section will 
confirm this result. 


2.4 Solution of the Auxiliary Integral Equation 

In the preceeding section, we stated that in the limit as the outer 
guide half-width became infinite, that (2.68) would become an integral 
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equation for the unknown expansion coefficient g fc . To this end, we 
rewrite the partial fraction expansion 


G(u) = 


+ I 


-'m 


to-8 , to+£ 

p m=l m 


(2.71) 


and the auxiliary matrix equation 


g = X 
s s 


+ I 




K -3 , 5 +5 

s p m=l s m 


( 2 . 68 ) 


As the outer guide recedes to infinite, we would expect 


G(to) 


to-B 


g.(, z I 

(0+Z 


dz 


(2.72) 


and 


g(w) = X (to) 


to— 3 


+ 


,si z l 

<o+z 


dz 


(2.73) 


A combination of (2.72) and (2.73) yields 


g(w) = A (to) G(<o) 


(2.74a) 


and 


G(w) = 



X (z) G(z) 

tO+Z 


dz 


(2.74b) 


Equation (2.74b) is an integral equation for G(to) when to is on the path 
o^o It can be shown that this expression is the auxiliary integral 
equation referred to in Section 2,3.2. 
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The expression for G(o>), given by (2.74b), has a branch singularity 
along the reflection of the path of integration. We evaluate G(oo) on 
each side of this branch cut and obtain 


G(-co| ) = -S’- + jir A (to) G(a>| ) + P.V, 

°3 P p °1 


M z ) G(0 . 

z-a) 


dz (2.75a) 


and 

iSsi. . (2.75b) 

Z-0) 

Subtracting (2,75b) from (2,75a), gives an expression for the discontinuity 
in G(cj) across the branch singularity 


G(— to _ ) = 


(o+S 


- JIT X(oj) G (to I ) + P.V. 


G (—to | ) -G (—co | ) = 2irj A (to) G(w| ) . (2.76a) 

°3 4 1 

This process for obtaining the branch discontinuity for a singular 
integral is based on the use of the Plemelj formula (Muskheliskvili , 
1964). A comparison of (2, 76a) with (2.48) shows that if G(w) is 
identified with T(m), that (2.74b) is the required auxiliary integral 
equation 


T(u>) - 


to— 0 


X (z) T (z) 
w+z 


(2.76b) 


In the following section, (2.76b) will be transformed into a suitable 
form for numerical methods of solution. 
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2,4.1 Conversion to a Smooth Kernel 

Close inspection of the partial kernel A(z) in equation (2.76b), 

-1 

reveals that the path passes dipectly over the poles of R (z). It 

will be shown in Section 2.5 that there is at least one zero of R(z) 

1/2 

lying in the open interval (k Q , k k Q ), w ^ en < > 1 and real. These 
zeros of R(z) correspond to TM-type surface modes within the dielectric 
slab. An attempt to solve for T (w) on the path will result in an 
integral of the Cauchy principal value type, which is in general, 
unsuitable for numerical methods of solution. 

Equation (2.761}), as it presently stands, has a singular kernel. 
This equation is converted to one which has a "smooth" kernel by 
changing the path of integration to one on which the kernel is analytic. 
Such a path is shown in Figure 6 and is labeled by a^. Equation (2.76b) 
then becomes 


T <“> -ir + ! dz (2 - 77) 

p ”1 

In continuously deforming the path from to a^, no singularities of 
the integrand of . (2. 76b) are encountered, since the integrand is analytic 
in region S(cf, Figure 6), Using the Cauchy theorem, the integrals in 
(2.76b) and (2.77) give identical results for to on the path 

Equation (2 .77) is no longer applicable when w is in region V of 
Figure 6, The effect of changing the path of integration from o ^ to 
is to shift the cut in the Rieman sheet for T(w) from a position 
directly above the path a ^ to one directly above a ^ (cf. Figure 6), 
However, we may analytically continue the solution for T(w) as given by 
(2.77) into region V by employing the relationship which describes the 
branch discontinuity, viz,, equation (2.48). The correct result for 




Figure 6. The integration paths for T(<d) and the region V of analytic 
continuation in the complex u> or z-plane. 
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T(oj), when u) is in region V, is given by 


T(o)) = 2irj X (— to) 


-1 
art - $ 


M g ) dz 

-w+z 



X(z) T (z) 
a)+z 


dz 


(2.78) 


Equations (2.77) and (2.78) are in a convenient form for numerical 
computation. We first solve (2.77) for T(to) on the path a^, and then 
substitute the proper values of oj on and a ^ into (2.77) and (2.78) 
respectively, to obtain T (to | ^ ) and T (—to | ^ ). 


2.4.2 Satisfaction of the Edge Condition 

In order that the solution obtained in Section 2.3 be acceptable, 
it must be demonstrated that the correct edge condition (Meixner, 1954) 
is satisfied at the guide edges, viz., (x = + b, z=0) . Due to the manner 
in which F(w) was factored in (2.35), it must be shown that T(w) 

Oi*r ; ] as I cu| -*■ «, where v is given by (2,39). For large |coj, (2.76) can 
be written in the following asymptotic form: 


T (co ) + 


MuiMai dz + 

u>+z 


(2.79) 


where A(z) is the asymptotic value of X(z), and is given by 


( 7 [Si] sin2(zb) s 1 - °) 




>. (2.80) 


0 [ exp (-z£ ) ] 


; £ > 0 


) 
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The limiting behavior of the infinite product expansions n( + uj,8) in 
(2.39) is deduced from an infinite product expansion of the gamma function 
(Whitehead, 1951). To check whether T(oj) has the correct asymptotic 
behavior, Z~ V ,,is substituted for T(z) in (2.79). The result is given 
by 



(2.81) 

For 0 < Re(v) < 1, these integrals can be evaluated by comparison with 
the appropriate Stieltjes* transforms (Bateman, 1954), and the results 
are given by 


OO 



-V 

7T CO 

sin(viT) 


(2.82) 



Substituting these results into (2.81), together with the assumed 
asymptotic expression for T(m), yields 
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-v * 

os -*■ < 


HI a 

_K+1J Sin(VTT) 


00 • 


St, = 0 


, l > 0 


(2.84) 


Comparing coefficients of w V in (2.84) gives the correct result (2.39), 
which is repeated here for clarity. 

\ 


1 - max 0 


, — arcsin - 

TT U 


K-l" 

<+l. 


/ 


; Z = D 


v = < 


> . 


(2.39) 


; SL > 0 


Thus, the correct edge condition is satisfied by the solution given in 
Section 2.3. 

2=5 Surface Waves in the Dielectric Slab 

The . existence of surface modes in the dielectric slab was previously 
mentioned as being attributed to the poles of R "^(a). In order to show 
the behavior of R(a) , we factor it into an equivalent form 

(2.85) 

( 2 . 86 ) 


R(a) = 7^ [(CO)-(CE)] 


where 


(CE) = t sinh(tt) + k£ cosh(xt) 


and 


(CO) = t cosh(xt)+ sinh(tt) . (2.87) 
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Equations (2.86) and (2.87) are respectively the determinantal or 
characteristic equations for the longitudinal propagation constants 
a for even and odd TM-type surface modes within the dielectric slab 
(Collin, 1960). 

We now write equation (2.1), which gives the magnetic intensity 
inside the slab, in terms of the even and odd characteristic equations 


(H ). 
y A 


1 

2 


k F (5 | a ) exp(-jab) 1 

it exp (50 
a 


cosh[t (z-A-t).] sinhft (z-ft-t) ] 

(CE) (CO) 

The integrand in (2„88) has simple poles at those values of a where (2.86) 
and (2.87) are zero and branch points at a = + k Q . The corresponding 
branch cuts of 5(a) are shown in Figure 7. The magnetic intensity of the 
surface modes can be evaluated by closing the contour of integration in 
(2.88) in the lower half a-plane for x > 0 and in the upper half a-plane 
for x < 0. Such closed contours are shown in Figure 7. It can be seen 
that in both cases; i.e., x > 0 and x < 0, the appropriate contour 
encloses all of the poles of the function R ''"(a). An application of 
the Cauchy theorem to the closed contour integrations yields the 
following results: (1) The radiation field can be identified with the 

contour integration along the branch singularity of 5(a). (2) The 

surface modes are identified with the residues of the integrand of (2.88) 
evaluated at the singularities of R ^(a). 

Before writing explicit expressions for the surface modes, some 
nomenclature for identifying each mode is presented. When k > 1 and real, 


j exp(-jax)da . (2.88) 
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equations (2.86) and (2.87) have a number of real zeros which lie in the 
1/2 

interval (k , k k ). We assign to the zeros of the even characteristic 
o o 

equation (2.86) even , integers , so that if there are Mkl such zeros, we 
obtain (oIq < < < a M) * For equation (2.87), the zeros are 

given by (ct^ < < . . . < %) • Moreover, the zeros of the even and odd 

characteristic equations are interleaved (Collin 1960) , and cIq is always 
the smallest of the set of values. We identify each surface mode with 
the corresponding value of.ou, i.e., the j-th mode is of the TMj o -type 
with respect to the x-directipn. If j is even or odd, we mean respectively, 
even or odd surface waves. Finally, the TMq surface mode. has no cut-roff, 
frequency, so that when k > 1, we always have at least one surface mode 
launched within the slab. 

We are now in a position to write the expressions for the magnetic 
intensity of the surface modes. When (m = 0,1,..., M) we express the even 
modes as 



For (n = l,3,o.o,N) the, odd modes are given by 



We see that both types of modes represent waves which travel away from 
the waveguide aperture in the x-direction and without attenuation. 
Furthermore, as the dielectric slab is moved away from the .waveguide 
aperture, the launching of surface waves becomes increasingly difficult, 
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This difficulty is due to the factor exp(££) in the denominator of (2.89) 
and (2.90), and represents an exponential damping of the surface mode 
excitation with increasing l. 

2.6 Radiation field jn the Open-Region 

The radiation fields in the open-regions C and E of Figure 2 may 
be evaluated from the integrations indicated by the corresponding 
Fourier transforms, viz., equations (2.3) and (2.5). However, the 
saddle-point method of integration (Collin, 1960) may be employed instead, 
and the radiation fields are obtained directly from the spectral weight 
coefficients C(y) and E(a). These weight coefficients need only be known 
in the visible range, i.e., 0 <_|a , yj<_ k Q . This limited knowledge is 
desirable from the standpoint of numerical computation. 

Before presenting the results of the saddle-point method of 
integration, we digress for the moment to investigate the effect that a 
dielectric or plasma slab has on a uniform plane wave. The result of 
this investigation will prove most useful in interpreting the equations 
for the radiated fields, 

2.6.1 Effect of the Slab on Plane Waves 

Before considering the expressions for the radiation fields, we first 
determine the influence that a dielectric of plasma slab has on a uniform 
plane wave. Figure 8 shows the incident field denoted by (E^, H^). The 
direction of propagation with respect to the normal of the slab is given 
by 9. The reflected and transmitted waves are.(E r , H^) and (E t , H t ) 
respectively, and the associated coefficients of reflectipn and 
transmission are given by R and T respectively. Also shown ip Figure 8 
are the forward and backward traveling waves (E^, H^) and E^» 
the slab, and the angle of refraction ip „ 


H^) inside 
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A lengthy but straightforward solution for the coefficients of . 
reflection and transmission yields 


and 


«“> - S& 

?(co = [RCcor 1 , 


(2.91) 


(2.92) 


where a = k Q sin(0). The functions, R(a) and Q(a), are identical 
to those employed in the preceding analysis and are given by (2.23) 
and (2.24) respectively. An inspection of the function Q(a) indicates 
that there are several angles of indidence 6 for which there is no 
reflection from the slab. Expanding (2.24) in terms of its zeros yields 


Q[k Q sin(9)] 


(k-1) (k+1) sin J^2t k Q »4:-sin^(0)J 

'' 3/2 ' — 

2j K cos(0) cosO);) 


[ sin(6) + [ sin(e) - J =5] . (2.93) 

The angles for zero reflection occur when k = l (slab removed) , at the 
Brewster angles 9_, (Jordan, 195Q) 

J3 


= + arcsin 



(2.94) 


and at certain other angles determined from the relation 


2t k 

o 


A 


K-sin (9) 


- 0 


sin 


(2.95) 
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The critical angles given by (2.95) are dependent not only on the 
relative dielectric constant of the slab, but also depend on the thickness 
An investigation of these angles shows that they correspond to a 
matched condition of the directional wave impedances as observed bn 
each side of the left-hand air-slab interface (cf. Figure 8). A simple 
example occurs when the angle of incident 0 is zero, and the slab 
thickness 2t is equal to an integral number of half-wavqlength as 

determined inside the slab. 

\ 

When the slab is a plasma, i.e., k < 1, we find no angles 0 where 
zero reflection occurs. Rather, the coefficient of reflection approaches 
unity as k + - °°, indicating transmission through the slab is cut-off 
for all angles of incidence. An asymptotic expression for the reflection 
coefficient is given by 


lim |R j 


|k 1/2 | cos(e) 

__ — — — — . — _> i 

|[2j 4- |< 1/2 (cos(0)]| 


(2.96a) 


On the other hand, the transmission factor T introduces A severe 
attenuation of the transmitted wave when the slab thickness is large 
or the relative dielectric constant is made negative. An asymptotic 
expansion for the coefficient of transmission is given by 

2 exp(-2t k o |ic 1/2 h 
| [2j + I* 1 ' 2 | cos(0)] | 


iim , Y | 

K->— 00 


(2.96b) 


Q[exp(-t |k 1/2 | )] 



We now have expressions for the coefficients of reflection and 
transmission for a uniform plane wave incident onto the slab. This 
information will be employed for a geometric interpretation of the 
solution of the problem at hand in the next section. 


2.6.2 Pattern in the Forward Direction 

The radiation field in region E of Figure 2 is obtained by applying 
a transformation of variable to the Fourier transform expression for the 
magnetic intensity. This transformation is given by 


I £ = j k Q cos(<}>) x = r sin(0) I 

a = k Q sin(<)>) (z - £ - 2t) = r cos(0) 


(2.97) 


Employing (2„97) in (2„5) yields 


( Ve 



exp[-j k Q r cos(0-<f>)] k Q Cos($) d<J). (2.98) 


Applying the method of saddle-point integration to (2.98) yields an 
asymptotic expression for the radiated field 


Urn (H ) 
r-x*> J 



cos 


(0) E[k o sin(0)J 


exp[-j (k^r-tr/4)] . 


(2.99) 


The radiation pattern in the forward direction, is obtained by 
substituting the solution for E given by (2.44) into (2.99) 
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PATTERN (0) « 



l[k sin(0) j , 


( 2 . 100 ) 


where F is the function constructed by the modified function-theoretic 
technique, viz., equation (2.35), and T is tdie trahsmission coefficient 
of the slab given by (2.92). The total power radiated in the forward 
direction is obtained from the expression 


P 


rad 


1 

(1+6°) b ir k 
p o 



F[jk Q cos(6)] 

R[k sxn(0)] 
o 


2 

de , 


( 2 . 101 ) 


where integration is only over positive angles of 0 due to the symmetry 
of the fields. The observation angle 0 and polar coordinates (r,9) for 
region E are shown in Figure 9. 

An examination of the pattern function (2.100) reveals that there 
are nulls in the forward pattern for those directions where F Jjk Q cos(0)J = 0. 
The angles 0 corresponding to these nulls can be identified with the 
directions of propagation in composite plane-wave rtepresentation of the 
waveguide modes (cf. Figure 10*). However, there is a maximum in the 
pattern at the angle corresponding to the incident mode exciting the guide. 

The salient features of the forward pattern can be summarized as follows: 


PATTERN NULLS: [0 - + arcsin(~-) ; s ^ p , s <_ N] (2.102) 

PATTERN MAX: [0 = + arcsinj^j ] , (2.103) 


where only propagating modes are considered, i.e., N <_ b k / ir. 
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2.6.3 Pattern in the Backward Direction 

The radiation pattern in the backward direction, region G, is 
obtained in exactly the same manner as was presented in Section 2.6.2 
for the forward region. The results are given by 

litn 
r-*°o 

exp[-j(k Q r - k Q b + tt/4)] (2.104) 



2 r cos (0 ) C[k Q sin(0)] ‘ 


PATTERN (0) = | F [jjk Q cos(0)J + 

R jk Q sin(0)J Fjjk Q cos ( Q )J exp £-2 j k Q ii. cos(0)JJ (2.105) 


and 


P 


rad 


1 

(1+6°) b it k 
P ° 



+ 


Q’ [k sin(6)] r -| I 2 

R- tk o 'Tin(0TT F [ jk o COs(e ^J I d& ‘ (2 ' 106) 

The angle of observation 0 for the backward direction is shown in Figure 
11. The pattern function, as given by (2,105), has a geometrical 
interpretation. We note that the radiation field is composed of two 
parts. The first is that due to F[-jk Q cos(0)], which we identify with 
the direct . field radiated from the guide aperture. The second part is 
due to F[jk Q cos(0)], which is the field radiated into the forward 
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direction. This forward field is modified by the function R, the 
reflection coefficient from the slab, and a phase difference 
exp[-2jk Q £ cos (6)], which is the path difference between the direct and 
reflected! fields. Figure 12 shows both fields being emitted from a 
common point, viz., the center of the waveguide aperture, and the 
path difference between the direct and indirect fields, i.e., u + v = 25,' 
cos(0) (cf. Figure 12 ). 
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3. SOLUTIONS FOR RELATED PROBLEMS 

In order to show the versatility of the modified function-theoretic 
technique, we present the solutions for two other open- region problems. 
These problems again represent a modification to the canonical problem, 
and their geometry is constructed by introducing a longitudinal 
inhomogeneity in one of the open- regions. 

The first of the modified problems is shown in Figure 13, and repre- 
sents an open-ended waveguide radiating into a reflecting plane. Such a 
geometry is closely related to large aperture reflector antennas. The 
solution for this problem is constructed from the function 

H(oj) = F 1 (oi) S(u) , (3.1) 

where F^ (to) is again the canonical function given by (2.36). The portion 
of H ( to) which represents the departure from the canonical function is 
represented by S(w). We require S(to) to have a simple 1 pole at w = 3^ 
and a branch cut along the path (cf» Figure 4). These requirements 
are identical to those for T(w), the modifying function for the problem 
in section 2. In fact, whenever there is a longitudinal inhomogeneity 
facing or opposite the waveguide aperture, the modifying function will 
have these requirements , 

A solution for the mode and weight coefficients in regions B and C 
are given by 

(-l) m H(-B I ) 
m'cr ' 

B — 

m b g (1+6°) 
s s 


(3.2) 
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PERFECT 

CONDUCTOR 


Figure 13. Waveguide radiating into a reflecting wall. 



Figure 14. The composite field for a waveguide radiating into a reflecting 
wall, u + v = 2i cos(e). 
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As for the problem in section 2, we see that the radiated field in 
region c is composed of two fields (cf. Figure 14). Referring to equation 
(3.4), we see that the direct field is given by H(-ui) . the indirect field 
is obtained by reflecting H(uO from the reflecting wall with its phase 
retarded by the path length y + v = 2 £cos(0). 

The geometry of the second problem we wish to discuss is shown in 
Figure 15. This problem represents radiation from a flanged waveguide. 
Since the longitudinal inhomogeneity is placed in back of the waveguide 
aperture -plane, we would not expect the modifying function to have a 
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branch singularity along the path a ^ (cf. Figure 4), and this is indeed 
the case. 

We again write the constructed function in the following form: 

G(u>) = F^aO U(u) : (3,8) 


We find that an analysis of the problem requires the modifying function 
U(o)) again have a simple pole at tu = however, for this geometry, 
the branch singularity is along the path instead of The solution 

for this problem is given as follows: 


- C-D s g<-6 8 | 03 > 

Bs b 6 (1+6°) 

s s 


(3.9) 


and 


E(0 = 


-jab r 


TT £ 


GUL ) + G(-Cl )e" 25d 


(3.10a) 


1 3 

The normalization condition is determined from the relationship 

G(B | ) = b g (-1) P (1+6°) . 

p'o p p 


(3.10b) 


The radiation pattern is given by 

PATTERN (0) = j G ^jk Q cos(0)J + 

G £-jk Q cos(0)J exp jj-2jk Q d cos(0)J | . (3,11) 

As in the previous examples, the radiation field can be interpreted in 

terms of the direct and indirect fields which originate from a common 

point (cf. Figure 16) . The expressions for the branch discontinuity and 

the auxiliary integral equation are as follows: 

U(uL ) - U(J ) = 2irj A(-u) U(-cd| ) 

1 °2 °3 


(3.12) 
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ahd 


U( W ) = 


f A(z) U(z) 


w-g 


Z+U) 


dz 


(3=13) 


The partial kernel \(%) is given by 
.jab 


A (-5) = 


a > 

sin(qb) exp(-2gd) °3 




(3.14) 


Another interesting problem is that of a flanged waveguide combined 
with a dielectric or plasma slab covering the guide aperture. Such is 
the case for radiation from an antehna flush mounted on the surface of a 
rocket with a plasma sheath covering the aperture-plane. The solution 
for this problem can be obtained by the modified function-theoretic 
technique. However, since the canonical problem is modified on both 
sides of the aperture-plane, we obtain two auxiliary equations. The 
form of the constructed function is given by 

J(co) = F 1 (w) [L(u)) + M(u>) ] . (3.15) 

The modifying functions, L(oo) and M(u) , have branch singularities along 
and respectively (cf. Figure 4). The simultaneous solution of 
the two auxiliary integral equations is not a simple one, since there 
is a coupling between L(w) and M(co) in the two equations. Hue to the 
complexity of this problem, we give no results, but only indicate that 
a solution is possible. 

A final note relating to the modified problems just discussed is 
that in each: case the correct edge condition may be shown to be 
satisfied. This is accomplished by comparing the asymptotic form of the 
auxiliary integral' equation with the appropriate Stieltjes transform 
(cf. Section 2.4.2). 



4. EXTENSION TO CIRCULAR GEOMETRY 
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The results that have been given thus far are concerned with the 
excitation of inhomogeneously filled open-regions by parallel-plate 
waveguides. Such problems are physically not realizable, but their 
solutions give a tremendous insight into the form of the solution for 
physically obtainable problems. Each of the problems discussed in the 
previous - sections can be solved by the modified function-theoretic 
technique. when the geometry is circular, and the electromagnetic fields 
have rotational symmetry. The only difficulty is that the factor N(w) 
employed in the construction of the canonical function F^(m), no longer 
has a closed form expression. We may, however, employ an integral 
formulation for N(co) given by Mittra and Bates (1965) 


N(oj) 


exp 


bo) 

IT 

i 



+ 


. bm 
J 2 


N 1 (m) 


(4,1) 


where 


lim 

N..(w) = a, c-*» 
1 a- c=b 


n(co,ri> 

n(u),s) 



i 



■^rj [7Tb lN ^ Czb) + J o (zb )) _1 i ^ 

° <4. 3) 


In Figure 17, we show the radii of the small and large waveguides, 
viz., b and a, for the associated closed-region problem, as the outer 
guide radius becomes infinitely large, (4.3) is determined as the limit 
of the quotient of the infinite product expansions for regions C and E. 
The functions J q and N q employed in (4.3) are respectively zero order 
Bessel functions the first and second kind. The function N^(w) is seen 
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to have a branch point at w = jk Q . The branch cut for N^(w) is chosen 

along the path (cf, Figure 18). 

Although (4.3) can not be written in a closed-form, we may 

approximate.it by a numerically truncated integral. We note that for 

large values of z the integrand of (4.3) is nearly equal to the 

corresponding. intergrand in the parallel-plate case, which does have a 

2 2-1 

closed-form expression. That is to say, [2/(iTbz) ] [N Q (zb) + J Q (zb)] -> 1 
for even moderately large values of z, e.g., this function equals 0.95 
when zb = 1. We, therefore, write (4.3) in the following approximate 
form: 


-b 


N 1 (u) s N 1 (w) exp < 2 ^ 


In 


1 - 




o 


:]K 


(zb) + N Q (zb) 


] 


-1 


- 1 


dz 


(4.4) 


Lit b z 

The function N^(oj) is the parallel-plate canonical function given by 
(2.38). The truncation point L is determined from a given numerical 
accuracy criterion. 

Aside from the integral expression (4.3) for N^(w), we are able to 
obtain two other useful relationships which aid in the numerical computation 
of this function. The method employed in obtaining these relationships 
is to find the limiting behavior of the product [n^(w)n^(-u>) ] as the 
outer waveguide radius becomes infinite. Here we have written n^ as 


the closed-region form of N.. 
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n t M - IW 
i n(w,5) 


n 

s=l* 

l 1 #A\ 

| exp | 

ui 

! W 

s o 

CO 

II a 8 
(-■ 

L u ^ 

exp 

( to 

1 / i ,2 

i /a -k 

l 3 ° 1 

, /i 2 -k 2 

[ s o 


(4.5) 


the values of y and a employed in (4.5) a ! re given by the zeros of the 
s s 

characteristic equations for regions C and E respectively (cf. Figure 17) 


J (y b) = 0 
o s 


J (a b) N (a a) - N (a b) J (a a) 
os os os os 


0 


(4.6) 

(4.7) 


Before taking the above indicated limit, we give two standard 
reflection formulas 


n(a3,n) n(-a),n) = 


i. [-ten 

« r ik \ z~\ 

3=1 L Vsl J 


J (db) N (a a) - N (ab) J (a a) 
o o o o 

J (k b) N (k_ a) - N (k b) J (k a) 
oo oo oo oo 


(4.8) 


and 


n(co f .O n(-« f O = 


oo 



2 1 

n 

1 - 

a 


S=1 

- 

Cfc 

S. 

_ 

oo 

““ 

fk \ 


n 

1 - 

0 

2 ] 

S=1 


a ‘ , 

8 / 



J Q (aa) 

J (k a) 
o o ' 


/ 2 2 

where a = /a) +k Q . Substituting (4.8), and (4.9) into the product 
[n^(m) n^(-w) ] yields 


J (ab) 
o 


[n 1 (w) n 1 (-w) ] = 


'N (aa) 
o 

J (aa) 

L. O 


- N (ab) 
o 


*7 (k b) 
o o 


r N (k a)" 
o o 

J (k a) 
L o o . 


- N (k b) 
o o 


(4.9) 


(4.10) 
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In the reflection relationships (4,8) and (4,9), we have employed the 
standard infinite product expansions for the appropriate Bessel functions 
(Abramowitz and Stegen, 1964) . Our aim is to find the relationship 
corresponding to (4.10) for the open-region case. To this end, we let 
a ->■ °° in (4.10) and obtain 


N (1) (ctb) 

N (co) N, (-(d) = - , (4.11) 

(kb) 

o o 

where and are Hankel functions of order zero, 

o o 

Equation (4.11) is quite useful in numerical calculations. The 

function must be known both on paths cr ^ and its reflection a^, shown 

in Figure 18. With the aid of this reflection formula (4.11), we need 

only numerically calculate N(co| ) from the integral (4.4), and obtain 

a 3 

N ( u) ] ) directly from (4.11). This process yields a fifty-percent 

1 

savings in integration time. Also, it is to our advantage to refrain 

from evaluating N( 0J | ) directly from the integral (4.4), due to the 

' a l 

logarithmic singularity of the integrand when w is on the path 

A second relationship involving the function N-^(u)) relates the 

function values on each side of the branch singularity. Such a 

relationship must be known when employing the function-theoretic 

technique. Evaluating (4.11) for on each side of the branch singularity, 

and by eliminating the common factor N(-w| ) , we obtain 

a 3 


N lHc 2 ) = 'Mo* 


r (ab) 

o 

N (1) (ab) 
o 


(4.12) 


where the paths a ^ and a 2 are shown in Figure 18. For large values of 
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to, (4,12) has the asymptotic representation 


lim 

oj-x® 


“N Coo| )- 
X 2 

V^r ) 

- i i . 


/ 2 ., 2 

/03 


= e 


-j (2 /a) +k b + tt/ 2) 
o 


(4.13) 


The corresponding relationship for parallel-plate geometry is given by 


N l (u) la 2 ) -j2 /j 2 +k 2 b 

T ) = 6 . (4.14) 

1 C 1 

Equations (4.13) and (4.14) show the close correspondence between the 
canonical function for circular geometry and for parallel-plane 
geometry in the asymptotic sense. 
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5 . NUMERICAL RESULTS FOR PARALLEL-PLATE GUIDE WITH SLAB 


Early attempts to solve the auxiliary integral equation (2.77) for 

T(oj| ,) were carried out by employing an iterative technique. The term 
°1 _ ± 

(z|^, - 3p) was substituted as a first approximation to T(z) in (2.77). 

The resulting numerical integration yielded a second order approximation. 
Subsequently, an n-th approximation was obtained by iteration of the 
(n-l)-th approximation. In general, it was found that for values of b 
and t which were of the same order of magnitude as the free-space wavelength 
Aq, the iterative method yielded a rapidly converging sequence of functions, 
i.e., only three or four iterations were required to give a solution which 
differed by less than one percent from the preceeding iterated result. 
However, when b and t were large compared to Aq, rapid divergence of the 
iterated result was encountered. Ih both the convergent and divergent 
cases, a plot of the logarithm of the magnitude of the iterated function 
for any fixed point on the path against the number of the iteration, 
yielded a straight line. Even the logarithm of the function values for 
the first and second iterations lay on that line, indicating that the 
spectrum of eigenvalues of the truncated auxiliary integral equation was 
of a special character. The character is that there exists one eigenvalue, 
dominant in modulus over all the others in the spectrum. The dominant 
character of this eigenvalue may be used to advantage in the construction 
of a solution for the divergent case. However, for divergent and slowly 
convergent cases, matrix inversion was employed instead. It was found 
that the computation time required to invert a 50 x 50 complex matrix was 


comparable to ten iterations. 
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Some numerical solutions to the problem shown in Figure 1 were 
obtained using a high speed digital computer, viz,, the IBM 360/75. In 
order to illustrate the effect of varying each of the physical parameters, 
the results -are divided into the following categories: (1) variation in 

b, the waveguide half-width, (2) variation in £, the slab -aperture 
separation, (3) variation in 2t, the slab thickness, and (4) variation 
in k, the relative permittivity of the slab. 

It was indicated in Section 2.6.2 that the number of nulls in the 
forward pattern is dependent on the waveguide half-width or in other 
words, equivalently dependent on the number of propagating guide-modes. 
Figures 19a through 19d show the radiation patterns for four different 
values of half-width. The values of b were chosen so that ail additional 
mode propagates in subsequent figure(s), beginning with one in 19a and four 
in 19d. In each ease, the incident mode was TEM, and this is the case for 
all results in this paper. Each pattern is normalized to unity in the 
broadside direction. Table 1 illustrates the distribution of power 
scattered into the various regions with the incident guide mode having' 
unit power. The various values of scattered power have been renormalized 
after computation in order that their sum is unity. However, the total 
scattered power before renormalization was found to be within four percent 
of unity in nearly all cases and occasionally not more than twelve percent. 
By increasing the number of sample points or the point of truncation in 
the numerical integration of the auxiliary integral equation, accuracy may 
be increased at the expense of computation time. However, it was found 
that such an increase in accuracy had very little effect on the "shape" 
of the radiation pattern, and only affected the distribution of scattered 
power. Typical machine-time for a single solution was about half-a-minute. 
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As indicated in Figure 19 and Table 1, the main lobe of the pattern 
narrows as the guide half-width is increased with little -af fect on the 
forward radiated power P £ . This result would also occur in the absence of 
the dielectric slab. An increase in b has a detrimental effect on the 
efficiency of surface mode excitation. Table 1 indicates that as the 
guide aperture is widened, the total power in the surface modes is 
decreased. 

The effect of changing the distance between the guide aperture and 
slab is next considered. Figures 20a through 20f show the radiation 
patterns when the distance t is varied from zero to A Q . As indicated in 
Section 2.6.2, the forward pattern is nearly independent of £, while its 
shape , depends on the transmission factor T. This independence of L is 
illustrated by the close similarity of the forward patterns of Figures 
20a through 20 f. The forward radiated power P^, shown in Table 2, is also 
nearly independent of t. 

The dependence of the pattern in the backward direction on t is 
not so simple, due to the phase relationship between the direct and re- 
flected fields in the composite field representation discussed in Section 
2.6.2, Table 2 compares the distribution of scattered power for various 
values of Z and shows that the excitation of surface modes becomes 
increasingly difficult as the slab is moved away from the waveguide 
aperture . 

In varying the slab thickness 2t, the pattern and scattered power 
have been obtained for various fractions of Aq, the free-spdce wavelength, 
instead of the wavelength within the slab. Figures 21a through 21f show 
the radiation patterns when the siab thickness varies from zero to A^. 




NJ > 


Table 2. The scattered power for various values of Z, the 
aperture-slab separation k = 1, k = 2, b = A 0 /4 
and 2t = A 0 /4. The free-space wavelength is A 0 . 
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Table 3, The scattered power for various values of 2t, the slab Width, 
given as fractions of A o , the free-space wavelength, 
k =1, k = 2, b = A /4, i, = A /4. 

O ’ o * o 
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watts 

watts 
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An interesting behavior of the total surface mode power is observed. As 
the slab thickness is increased, the amount of power carried by the surface 
modes increases. However, when the slab is sufficiently thick, allowing 
another mode to propagate, the total power decreases . A further increase 
in slab thickness again increases the surface mode power until another 
mode begins to propagate. 

Numerical results for various values of slab permittivity are shown 
in Figures 22a through 21e. radiation patterns are shown only for positive 
values of K, since little power is radiated when the slab is a plasma 
and when the aperture-slab separation £ is zero. The effect of an increase 
in k is to broaden the forward lobe of the pattern. This effect is 
accompanied by a decrease in the forward radiated power P , indicated in 

£j 

Table 4. Also shown in Table 4 are several cases for plasma slabs. When 
k < Oj a large fraction of the power is reflected back into the waveguide, 
as indicated in the column headed Pg. For large negative values of k, 
both the forward and backward radiated power is reduced; however, there is 
a greater fraction of power leaking around the guide edge into region C 
than is radiated through the slab into region E. This is the case when 
£ = 0. When the plasma slab is not directly against the guide aperture, 
an appreciable amount of power will be radiated into the backward region. 
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Table 4. The scattered power for various values of k, the relative 
dielectric constant of the slab, k * 1, b » A 0 /4, 2t = 
A 0 /4, l - 0. The free-space wavelength is A 0 . 
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6. CONCLUSIONS 

The present solution for radiation from an open-ended parallel- 
plate waveguide through a dielectric or plasma slab lends itself to a 
numerically accurate and rapid solution. Typically, the function-theoretic 
technique and its subsequent extensions and modifications (hereafter called 
MFT) have several distinct advantages over other methods, compared with 
some variational techniques for the same problem, the MFT yields at least 
the most useful information, viz., waveguide reflection coefficients. 
However, with little additional but straightforward computations, all of 
the scattered fields can be determined. Formulation in terms of the 
characteristic modes for each region in the MFT leads to a simple mode 
matching procedure. The process of continuing the incident waveguide 
field into the open-region, as is the case when the Weiner-Hopf technique 
is employed, is not in line with a physical understanding of the problem. 
Finally, the correct edge condition is automatically satisfied, since 
a priori information of the asymptotic behavior of the fields in the 
neighborhood of the waveguide edge is incorporated in the MFT. In 
contrast, the satisfaction of this important condition is difficult - if 
not impossible - to show when other methods are employed. For example, 
Mittra (1963) demonstrated that satisfaction of the edge condition in 
the solution of a doubly-infinite set of matrix equations relating to the 
closed-region waveguide bifurcation problem, was strongly dependent upon 
the way in which the number of equations was determined in the truncation 
procedure. 

Application of the MFT to waveguide problems other than phe basic 
parallel-plate geometry, requires only a straightforward extension of 
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the method employed for the simpler geometry. The solutions for 
circular guides with rotationally independent fields are easily obtained, 
and they possess all of the desirable characteristics of parallel-plate 
geometry. Also the MFT is easily extended to geometries containing layered 
media in the open— regions , Such layered problems are encountered in 
diagnostics of laboratory plasmas, where a bulk plasma is surrounded by 
an ion sheath and glass envelope. 

Other problems which are easily attacked by the MFT include: (1) 

The flanged waveguide; (2) radiation from an open-ended waveguide onto 
a reflecting screen; and (3) a flanged guide radiating through a plasma or 
dielectric slab. In fact, most open-ended waveguide problems in which 
the geometry is uniform in the transverse direction but longitudinally 
non-uniform can be formulated and solved by the MFT. 
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APPENDIX 


ESTIMATION OF THE TRUNCATION ERROR FOR INFINITE PRODUCT EXPANSIONS 


In the numerical solution of problems via the function-theoretic 
technique and its subsequent modifications, one is required to evaluate 
infinite product expansions of the following form: 


n(w,e) = n 

S=1 


exp'(Ve g ) 


where g ** v pr-\ ^ , k and b are constants described 

s \bj o*o 

elsewhere, and w is the function variable. To an end of evaluating (Al) 
in a truncated form, together with an estimate of the truncation error, 
we write 


JI(u),g) = I IT (1 - w/g g ) exp(oj/g g ) e R 
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The expression for the truncation error is given by 


= exp I [ln(l - cu/3 ) + u)/g ] 
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Expanding the logrithmic part of (A3) in a power series when 
(K + 1) ir/(2bk ) >1 u)/k I > 1, we obtain 
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Interchanging the indices of summation, i.e», summing on t first, together 


with the relationship 
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which is valid for (K+l) > y 


truncation error for rn-real 
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(A5) 


, we obtain bounds for the 
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The summations in (A6) can be evaluated exactly (Jolly, 1961) 

t 


I 

t=2 


x 


t(t-l) 


= x + (1-x) ln(l-x) 


(A6) 


(A7) 


Substituting (A7) into (A6) and taking only the first term in the power 
series expansion of In(l-x) , we find that e has the following asymptotic 
representation: 



t 


(A8) 


when K >> <, 

IT 

Equation (A8) is found numerically to yield excellent results for 
M < 10k o and K = 100c To a first order approximation, the value of 
e is the same for both + w. One may check the accuracy of (Al) as 
computed from (A2) by use of the following reflection formula in which 
the truncation error is cutaulative: 
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k sin(ab) 

n( w ,3) n(- u ,e) - ~d£aTb) '» (A9) 

n 2 

where a = /k +to . We further note that in the solution of the 
o 

auxiliary integral equation (cf. section 2.4), factors of the form 
n(to,6)/II(-co,3) are encountered. Due to the even symmetry of the truncation 
error (A8) , we see that in the above quotient of infinite product 
expansions, the truncation error is self-canceliiig. 
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